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The thermal conductivity k is measured in a series of Cui-^Mg^GeOa single crystals in magnetic 
fields up to 16 T. It has turned out that heat transport by spin excitations is coherent for lightly 
doped samples, in which the spin-Peierls (SP) transition exists, at temperatures well above the SP 
transition temperature Tsp. Depression of this spin heat transport appears below T* , at which the 
spin gap locally opens. T* is not modified with Mg-doping and that T* for each Mg-doped sample 
remains as high as Tsp of pure CuGeOs, in contrast to Tsp which is strongly suppressed with 
Mg-doping. The spin-gap opening enhances phonon part of the heat transport because of reduced 
scattering by the spin excitations, producing an unusual peak. This peak diminishes when the spin 
gap is suppressed both in magnetic fields and with the Mg-doping. 

PACS numbers: 66.70.+f, 75.30.Kz, 75.50.Ee 



I. INTRODUCTION 

The discovery of the first inorganic spin-Peierls (SP) 
compound CuGeC>3 (Ref. |]) has triggered subsequent 
studies of the, impurity-substitution effect on the spin- 
singlet states,H and the existence of the disorder-induced 
transition into three-dimensionajjfintifcrromagnctic (3D- 
AF) state has been establishedBIa Since the relevant ex- 
change energies, i.e., intra- and interchain coupling J and 
J', do not change with doping except at the impurity 
sites, it is difficult to understand the impurity-induced 
transition from the SP to 3D-AF state in the framework 
of the "conventional" competition between dimensionali- 
ties (where J'/ J is an essential parameters) and therefore 
such impurity-substitution effect is a matter of current 
interest. 

In CuGeC>3, a small amount of impurity leads to an 
exotic low-temperature phase where the lattice dimeriza- 
tion and antiferromagnetic staggered moments simulta- 
neously appear [dimerized AF (D-AF) phase] B Moreover, 
when the impurity concentration x exceeds a critical con- 
centration x c , thefSP transition measured by dc suscep- 
tibility disappearsQia and a uniform AF (U-AF) jphase 
appears below the Neel temperature Tn ~ 4 K.B The 
D-AF ground state can be understood as a state of spa- 
tially modulated staggered moments accompanied with 
the lattice distortion.EE3 However, the mechanism of the 
depression of the SP phase and the establishment of the 
disorder-induced antiferromagnetism are still to be elu- 
citated. A transport measurement could be a desirable 
tool in dealing with such a problem, because the mobility 
of spin excitations or spin diffusivity is often sensitive to 
the impurities. Nevertheless, even a crude estimation of 
the mobility has not been carried out for the spin exci- 
tations so far. 

Recently, we reported intriguing behaviors of the ther- 
mal conductivity k of pure CuGeOsJiil In Ref. [ll], the 



existence of the spin heat channel (k s ) was suggested. 
However, we cannot separate k s and the phonon thermal 
conductivity K p h with the measurement of the pure sam- 
ple only. In this work, we have resolved this problem with 
Mg-doped crystals. It turned out that the spin diffusion 
length, which can be calculated from k s , is much larger 
than the distance between adjacent spins for the pure 
CuGeC>3, indicating coherent heat transport due to the 
spin excitations. Since the spin heat transport almost 
disappears for the heavily doped samples in which the 
long-range SP ordering is absent, it is suggested that the 
large mobility of the spin excitations plays an important 
role in the SP ordering. 

The spin heat transport is also a probe of the spin gap. 
It is found that local spin-gap opening is robust against 
the impurity doping, though the temperature of the long- 
range ordering is strongly suppressed with x. Phonon 
heat channel (K p h) provides an unusual thermal conduc- 
tivity peak in the SP-state, which was one of the topics of 
our previous report S3 The peak is drastically suppressed 
with both magnetic field and Mg-doping, indicating that 
the spin-gap opening has produced the peak in K p h via 
strong phonon-spin interaction. 

II. EXPERIMENTAL 

The single crystals of Cui^^Mg^GeOa were grown with 
a floating-zone method. The Mg concentration is de- 
termined by inductively coupled plasma-atomic emission 
spectroscopy (ICP-AES).El The critical concentration x c 
was carefully determined by Masuda _ei al. for the same 
series of crystals that we have usedfl'tj They found that 
the Neel temperature Tjv jumps at the impurity-driven 
transition from D-AF to U-AF phase and that phase is 
separated in the transition range of 0.023 fjnp]< 0.027, 
indicating that the transition is of first order .BIL3 We used 
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six samples with x = 0, < x < x c , and x > x c for the 
thermal conductivity measurement. They were originally 
prepared for the neutron and the synchrotron x-ray ex- 
periments, and thus the sample homogeneity has been 
confirmed in detail which is described elsewhere.113 Typi- 
cal sample size is 0.2 x 2 x 7 mm 3 (a x b x c). Since our 
interest is thermal conductivity along the ID spin chain, 
we need the sample with a longer dimension along the 
c axis. Thus we payed particular attention to the Mg 
homogeneity along the c axis for this study, and the best 
way is to use the samples prepared for the neutron ex- 
periments. The fundamental magnetic properties-of the 
samples used already appeared in another paper .tn 

The thermal conductivity is measured using a "one 
heater, two .thermometers" method as the previous 
measurement O The base of the sample is anchored to 
a copper block held at desired temperatures. A strain- 
gauge heater is used to heat the sample. A matched 
pair of microchip Cernox thermometers are carefully cal- 
ibrated and then mounted on the sample. The maximum 
temperature difference between the sample and the block 
is 0.2 K. Typical temperature difference between the two 
thermometers is ~0.5 K above 20 K and 0.2 K below 
20 K. The temperature-controlled block is covered with 
a metal shield so that heat loss through radiation be- 
comes negligible. Magnetic field is applied parallel to the 
c-axis direction. The measured temperature range is be- 
low 30 K, which is well below the temperature scale of 
the intra-chain coupling J c 120 K according to neu- 
tron scattering measurement!! 3 !) , and therefore, the one- 
dimensional quantum spin liquid (1D-QSL) is realized 
above the transition temperatures. 



III. RESULTS 



Figure 1 shows the temperature dependence of k for 
samples with and without Mg substitution: l(a)-l(f) are 
for x = 0, 0.016, 0.0216, 0.0244, 0.032 and 0.040, re- 
spectively, x = 0.016 is well below x c , x = 0.0216 is 
close to x c , x — 0.0244 is on the border and x = 0.032 
and x ~ 0.040 are in the U-AF region. The transition 
temperatures T£p for the x < x c samples and for 
the x > x c samples, ssihich were determined by suscep- 
tibility measurement ,113 are indicated by arrows in the 
figures. A sharp drop is observed just below T§ p for the 
pure CuGeC>3 as presented in the inset of Fig. 1(a). The 
x = 0.016 sample shows the same behavior at Tsp, while 
n is enhanced just below Tgp for the x = 0.0216 and 
0.0244 samples. 
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FIG. 1. Thermal conductivity of Cui-^Mg^GeOa single 
crystals in magnetic fields up to 15 or 16 T applied along the 
chain direction, for (a) x — 0, (b) x = 0.016, (c) x = 0.0216, 
(d) x = 0.0244, (e) x = 0.032, and (f) x = 0.040. Inset: a 
magnified view of k near the spin-Peierls transition in zero 
field. The dashed line is an extrapolation of k from above the 
transition temperature 

In Figs. 1(e) and 1(f), a slight change in the curvature 
is observed at Tjy and the field dependence in the Neel 
state is opposite to that above Tjv ; K decreases with field 
below Tjv while k increases with field above Tjy. The 
three-dimensional magnon excitations, which appear in 
the Neel state, are responsible for this behavior. The 
thermal conductiidiy in U-AF phase will be further dis- 
cussed elsewhereEj Although the Neel transition (to D- 
AF phase) should be present also for the x < x c samples, 
we cannot identify the corresponding feature at the tran- 
sition, because T/v (~ 2 — 2.5 K)cl is too close to the lower 
temperature limit of our measurement, ~ 2 K. 

Taking a look at the full scale of each figure, one can 
notice that k is strongly suppressed with Mg substitution 
in the whole temperature range. Two thermal conductiv- 
ity peaks are present in the zero-field curve of Fig. 1(a); 
one is in the SP phase and the other is in the 1D-QSL 
state. The broad shape of the high-temperature peak 
is maintained with doping. On the other hand, the low- 
temperature peak (SP peak) diminishes with x < x c , and 
disappears for the x > x c samples. 

The SP peak is suppressed also in magnetic fields 
[Figs. l(a)-l(d)]. In order to see the field dependence 
in more detail, the thermal conductivity at the SP peak 
is plotted as a function of magnetic field in Fig. 2 for 
the x — sample, n is suppressed with magnetic field 
and suddenly drops at the threshold field, H c , at which 
the system makes a transition to the incommensurate 
phase. The same features were observed in our previous 
measurement with the crystal grown by the laser-heated 
pedestal growth technique .liil Also, the typical value of k, 
for example at the high-temperature peak, is almost the 
same as that of the previous result. We thus believe that 
the above features are intrinsic in this material regard- 
less of the growth method ofpthe crystal, as long as high 
quality single crystal is used.E-3 
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FIG. 2. Magnetic field dependence of the thermal conduc- 
tivity at 8 K for pure CuGeOs, demonstrating the suppression 
of the low-temperature peak with field application. 

In a temperature range close to Tj? p , K increases with 
field for the Mg-doped samples, while little field depen- 
dence is seen in the pure sample [Figs. l(a)-l(d)]. Note 
that the field dependence extends even above T^ p for 
the Mg-doped samples. The variation in the field depen- 
dence implies qualitative difference of the SP transition 
between pure CuGe0 3 and doped CuGe0 3 . This field 
dependence is apparent also in the x — 0.032 sample, 
which is above x c [Fig. 1(e)], but rapidly diminishes in 
the heavily doped sample [Fig. 1(f)]. 

IV. DISCUSSION 

CuGeOs is an insulator and heat conduction by elec- 
trons or holes is absent unlike metals. Instead, low- 
energy spin excitationala can carry heat. In the 1D-QSL 
state, total k is given by a sum of k s and « p h as, 

k = k s -\- K ph (in ID - QSL, T SP < T <C J). (1) 

In contrast, k s is suppressed at the lowest tempera- 
ture region in SP phase because few spin excitations are 
present due to the spin gap, and the total k represents a 
phonon contribution only, 

K~Kph (T«T SP ). (2) 

The field and impurity dependence of K p h below Tsp will 
be discussed in subsection A. The value of n s in 1D-QSL 
will be crudely estimated in subsection B. Finally, the 
problem of local spin-gap formation will be dealt with in 
subsection C, using this probe. 

A. k below Tsp: depression of the SP peak in 
magnetic fields and with Mg-doping 

Since the low-temperature peak is drastically sup- 
pressed with the application of field as shown in Figs. 1 
and 2, the feature should be related to the spin-excitation 
spectrum. Wp-have proposed the following explanation 
for this peak.Eil At temperatures well below Tgp, where 
low-energy spin excitations are negligible, heat is carried 
mostly by phonons and K p h increases with temperature 
because of the growing population of phonons like usual 
insulating crystals. With increasing temperature, the 
increasing number of thermally excited spin excitations 
scatter phonons, and K p h diminishes. As a result, the 
thermal conductivity peak shows a peak. Although the 
thermal spin excitations begin to carry heat, the increase 
in k s appears only in the temperature region slightly be- 
low Tsp [see inset in Fig. 1(a)]. 

In order to confirm the above understanding, the impu- 
rity substitution is helpful because the spin gap is absent 
in heavily substituted samples. We normalized the size 
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of the peak for each sample as k(5 K)/k(15 K), and plot- 
ted it against x in Fig. 3. One can see the suppression of 
the peak with x similarly to the case of field application 
(Fig. 2). Moreover, the peak disappears at the concen- 
tration x > x c where the long-range SP order no longer 
exists, as shown in Figs. 1(e) and 1(f). We have thus 
confirmed the correlation between spin gap and the peak 
also by Mg substitution, giving the additional evidence 
for the scenario described in the previous paragraph. 
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3. A normalized size of the low-temperature peak, 
k(5 K)/k(15 K), is plotted against Mg-concentration. A sup- 
pression of the peak with Mg doping is evident. 



The x = 0.0244 sample, which is on the border from 
SP (D-AF) to U-AF demonstrates enhancement of k be- 
low Tg P , as observed in Fig. 1(d). Since the enhancement 
is suppressed in magnetic fields, the existence of a well- 
defined spin gap is suggested even at x — x c . We can 
notice the rapid change in k s in the (low-temperature) 
vicinity of Tgp, which is another sign of the spin-gap 
opening, for the x = and 0.016 samples [Figs. 1(a) and 
1(b)]. However, this feature is absent for the x = 0.0216 
and 0.0244 samples [Figs. 1(c) and 1(d)], since k s is sig- 
nificantly reduced with x because of impurity scattering. 



Figure 4(a) shows the zero-field thermal conductivity 
of all the samples. The thermal conductivity is strongly 
suppressed with Mg-substitution in the whole measured 
temperature range. The x dependence of k below ~4K 
can be understood as the difference in the scattering rate 
of phonons. Since we showed in Ref. |ll| that the scatter- 
ing by planer defects is dominant for the x = sample, 
difference in the number of planer defects may cause the 
variation of K p h in our series of samples. The number of 
planer defects does not necessarily increase with increas- 
ing x; note that K in the x = 0.032 sample is smaller than 
that in the x = 0.040 sample below ~ 5 K. 
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FIG. 4. (a) Temperature dependence of the zero-field ther- 
mal conductivity of Cui-^Mg^GeOs single crystals, (b) The 
thermal conductivity at 17.5 K as a function of Mg concen- 
tration. The spin-Peierls temperatures measured by dc sus- 
ceptibility are plotted together. 



B. k above Tsp: 



coherent heat transport in the 
1D-QSL 



Remembering the origin of the low-T thermal conduc- 
tivity peak, we can notice that scattering by the spin 
excitations, in turn, becomes dominant in K p h above the 
peak temperature for the x < x c samples. Noting that all 
the x > 0.016 curves looks parallel to one another above 
~ 15 K, Kph in the x > x c samples should probably have a 
temperature dependence similar to that in the x = 0.016, 
0.0216 and 0.0244 samples above ~ 15 K. Therefore, it 
is expected that phonon heat transport is governed by 
the spin-scattering also for the x > x c samples, at least 
above ~ 15 K. 

In contrast to rather complicated temperature depen- 
dence below Tsp, k(T) in the 1D-QSL region is simpler, 
which is an advantage of discussing x dependence in the 
region. As a crude approximation, it is assumed that 
the phonon part is not so much x dependent and that 
most of the ir-dependenee. comes from k s , because the 
detailed characterizationo guarantees good homogene- 
ity in both pure and Mg-doped crystals and the direct 
modification in the phonon modes due to the Mg substi- 
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tution can be estimated to be negligible. The scattering 
rate of phonons due to the point defects is more than 
two orders of magnitude smaller than the total scatter- 
ing rateO One may notice that the robustness of K p h to 
the Mg-doping is natural, considering that phonons are 
mainly scattered by spin excitations above ~ 15 K, which 
is the conclusion of the previous paragraph. Since the 
population of both phonons and spin excitations does not 
change so much with x in the temperature rangepabave 
Tsp, as has been reported in specific heat results,Qii30 x 
dependence of K p h is expected to be small in the 1D-QSL 
state.E9 

Figure 4(b) shows the x dependence of k at 17.5 K. k 
rapidly decreases with x up to x c , and it saturates when 
x exceeds x c . Therefore, it is naturally assumed that k s 
is close to zero in the heavily doped samples with x > 
x c . This result suggests some correlation between the 
mobility of the spin excitations and the SP ordering. 

Following the above interpretation, it is possible to es- 
timate k s by subtracting k of the x — 0.040 sample from 
the measured k for each sample. The result indicates that 
most part of the k of pure CuGe03 is due to n s , which is 
approximately 0.4 W/cm K at 17.5 K. This gives the dif- 
fusivity of the spin excitations-jDs of ~1.6xl0~ 3 m 2 /s, if 
we assume C s ~1.8 J/mol K.LLH D s in ID systems can be 
written as v s l s , where v s is the velocity of the spin excita- 
tions and l s is the mean free path, in other words, the spin 
diffusion length. Assuming the spittwave like dispersion 
given by des Cloizeaux and Pearsonll3 e=(n/2) J c \ sin(fcc)| 
and J c ~ 120 K,E^I one can evaluate v s ~ (n /2)cJ c /h= 
1.3xl0 4 m/s and l s ~ 130 nm (e is the energy of the 
spin excitations, k is the wave vector and c (~ 0.3 nm) 
is the distance between adjacent spins). Note that the 
spin diffusion length of the same order of magnitude was 
obtained from the NMR relaxation measurement for an- 
other ID spin system, AgVP2S6-til 

Since l s /c ~ 430 is much larger than unity, we can con- 
clude that the spin heat transport for the pure CuGeC>3 
is coherent. If the same estimation is applied also to 
the Mg-doped samples, l s /c ~ 100 for the x = 0.016 
(< x c ) and l s /c ~ 50 even for x = x c at this tempera- 
ture. Since c/x c is approximately calculated as 40, the 
above estimation indicates that the spin diffusion length 
in the 1D-QSL state exceeds the mean impurity distance 
c/x, as long as x is less than x c . Whereas, the spin exci- 
tations are not so mobile when x > x r . 



C. k just above T^ p : 



short-range SP order 



Figure 5 shows k — k(x — 0.040) (= Ak) of the three 
Mg-doped samples which have SP transition. One can 
find that An is almost temperature independent above 
T* ~ 15 K and that Ak of all the samples deviates from 
this temperature-independent value below T*. Since T* 
is close to Tsp of the pure CuGeOs, the deviation can 
be attributed to a precursor of the SP transition. Note 



that similar behavior is observed in the Raman scattering 
measurement on Zn- and Si-dopd CuGe03.t3 
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5. Temperature 
pendence of Ak [= k — k(x = 0.040)], for the x — 0.016, 
0.0216 and 0.0244 samples. The local spin gap opens below 
~ 15 K in all the samples. 



Recently, it is shown by synchrotron x-ray diffraction 
measurement that FWHM of the Bragg peak from the 
lattice dimerization reaches the resolution-Limit at a tem- 
perature Tgp~ ray which is far below TgpE3 The authors 
claimed that the long-range order takes place only below 
Tgp~ lay and— that lattice dimerization is short-range just 
below Tg^p.LJ The gap-like feature in k is also to be ex- 
plained along the idea of the short-range order (SRO). 
However, the data in Fig. 5 suggests that SRO grows 
from T* ~ 15 K, which is even above Tgp. 



The field dependence of k is a strong indication of the 
presence of the SRO below T* . Figure 6 shows the tem- 
perature dependence of Ak for Mg-doped samples in vari- 
ous magnetic fields. In all cases, we observed the recovery 
of Ak with increasing fields below T*. The results are 
consistent with the notion that the reduction of Ak below 
T* is due to the development of the local spin gap. Since 
the magnetic field is thought to reduce the magnitude of 
spin gap even above T,?p, the number of spin excitations 
responsible for the heat transport will increase with the 
field. Moreover, the field dependence only appears be- 
low T* in our observation. We have no field dependence 
above T*, which is consistent with the notion that even 
the local spin gap is absent. 
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FIG. 6. Temperature dependence of Ak for the (a) 
x = 0.016, (b) x = 0.0216, (c) x = 0.0244 and (d) x = 0.032 
samples. Magnetic fields are applied along the chain direc- 
tion. 



It should be emphasized that the heat transport by 
spin excitations is very sensitive to the spin-Pcicrls SRO. 
We observed a magnetic field dependence in Ak below 
T* even for the sample with x — 0.032 (> x c ), where the 
spin-Peierls LRO no longer exist in the whole tempera- 
ture range. According to the synchrotron x-ray measure- 
ment, SE-SRO develops below 7.5 K in the x — 0.032 
sample.E 3 ] For the x = 0.040 sample, however, there was 
no evidence that SP-SRO actually develops. In contrast, 
one can see that the magnetic field dependence of k still 
exists for the x — 0.040 sample (Fig. 1(f)) even though 
the change is not so obvious as in the x = 0.032 sam- 
ple. We think that the difference of a length scale of 
probed area between the two measurement techniques 
causes such discrepancy. 

The above results tell us that the impurity-substitution 
modifies the ground state and suppresses the spin-Peierls 
ordering in CuGeOa system through different mecha- 
nisms. The importance of the inter-chain coupling and 
the tendency to the 3D-antiferromagnetism increases be- 
low ~ Tjv. As a result, Tjv increases with ai_and the 
ground state changes from D-AF to U-AF at x c .Efl On the 
other hand, the local spin-gap formation, which occurs at 
T* , is goverened only by one-dimensional nature of the 
spin system coupled with 3D phonons and has nothing to 
do with 3D-AF fluctuation. Therefore, the temperature 
of SP-SRO is not modified with long as the spin- 

singlet state is energetically favored, (the spin-singlet 
state will be no more favored in a heavily disordered 
system where long- wave-length spin excitations, whoae 
energy is less than the spin-gap energy, are absent EB) 
The spin-diffusion length l s , in turn, rapidly diminishes 
with x owing to the growing impurity-scattering of the 
spin excitations. Since the spin correlation length is di- 
rectly related to l a ,o the length scale of the SP domains 
is reduced with x. As a result, signals of the spin gap, 
detected with any probe, diminishes in size with x. 



V. SUMMARY 

The thermal conductivity of the Mg-doped CuGeOa 
helps to understand both the spin-gapped state below 
the SP transiton and the 1D-QSL state above the tran- 
sition. Large spin heat transport is observed in the 1D- 
QSL state for pure CuGeOa, and rapidly diminishes with 
Mg-doping, accompanied by the suppression of the long- 
range SP ordering. The spin gap opening suppresses k s 
and enhances /t p h- Examining the x dependence of the 
spin-gap features, it turned out that the local spin gap 
opens at a temperature (T*) independent of x, suggest- 
ing that the suppression of SP ordering is attributed to 
the reduction of the spin diffusion due to the impurity 
scattering. It is expected that the above analysis of the 
impurity-substitution effect on the thermal conductivity 
is applicable to other one-dimensional spin systems and 
that thus obtained transport properties may reveal new 
aspects on such materials. 
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APPENDIX: PHONON SCATTERING RATE 
DUE TO POINT DEFECTS 

We can show that the scattering rate 1 /r p by point de- 
fects, which is introduced by the substitution with the ion 
with different mass, is very small even for the most heav- 
ily doped sample (aj— j= 0.040). The dominant-phonon 
approximation gives J13 



na 3 (k B T) 4 



4rwre 3 ?T(AM/M) s 



(Al) 



[n (= x = 0.040) is density ratio of point defects, a 
(~4x 10~ 8 cm) is lattice constant, m (= 15) is number 
of phonon modes (3 times the number of atoms per unit 
cell), v (~ 5 x 10 5 cm/s) is phonon velocity and AM/M 
(= 0.62) is the ratio (mass difference between Mg and Cu 
atoms) /(mass of Cu atom)]. l/r p can be estimated to be 
about 1.0 x 10 s (s^ 1 ) at 30 K. 

In comparison, the total scattering rate 1/r for the 
x = 0.040 sample can be obtained from the data as 



1/r 



a 



3ftph 



(A2) 



assuming phonon specific heat C p h to be ~ f3T 3 and us- 
ing the value of [3 ~ 2.8 x 10~ 6 (J/cm 3 K 4 ), shown in 
Rcf. |ll|. Taking K p h ~ 0.1 (W/cm K) from our result of 
the k measurement, 1/r ~ 6.4 x 10 10 (s _1 ) at 30 K. Since 
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this 1 /r is more than 600 times larger than 1 /t p (the ratio 
is even larger below 30 K), the scattering by the point de- 
fects is not important for the three-dimensional phonon 
heat transport in the measured temperature range. 
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